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Abstract 

We consider a particle diffusing along the links of a general graph possessing some absorbing 
vertices. The particle, with a spatially-dependent diffusion constant D(x) is subjected to a 
drift U(x) that is defined in every point of each link. 

We establish the boundary conditions to be used at the vertices and we derive general 
expressions for the average time spent on a part of the graph before absorption and, also, 
for the Laplace transform of the joint law of the occupation times. Exit times distributions 
and splitting probabilities are also studied and several examples are discussed. 

1 Introduction 

For many years, graphs have interested physicists as well as mathematicians. For instance, 
equilibrium statistical physics widely uses model systems denned on lattices, the most popular 
being certainly the Ising model [lj. On another hand, in solid-state physics, tight-binding 
models (see, for instance, [2]) involve discretized versions of Schrodinger operators on graphs. 
For all those models, the thermodynamic properties of the system heavily depend on geometrical 
characteristics of the lattice such as the connectivity and the dimensionality of the embedding 
space. However, in general, they don't depend explicitly on the lengths of the edges. Random 
walks on graphs, where a particle hops from one vertex to one of its nearest-neighbours, have 
also been studied by considering discrete Laplacian operators on graphs [3]. 

Such Laplacian operators can also be useful if they are defined on each link of the graph. For 
example, in the context of organic molecules f4], they can describe free electrons on networks 
made of one-dimensional wires. Many other applications can be found in the physical literature. 
Let us simply cite the study of vibrational properties of fractal structures such as the Sierpin- 
ski gasket [5j or the investigation of quantum transport in mesoscopic physics 0[7]. Weakly 
disordered systems can also be studied in that context [8]. It appears that weak localization 
corrections in the presence of an eventual magnetic field are related to a spectral determinant on 
the graph. This last quantity is actually of central importance and interesting by itself [9| 110] . 
In particular, it allows to recover a trace formula that was first derived by Roth [TT]. Moreover, 
the spectral determinant, when computed with generalized boundary conditions at the vertices, 
is useful to enumerate constrained random walks on a general graph |12| . a problem that has 
been addressed many times in the mathematical literature |13j . 



Brownian motion on graphs is also worthwhile to be investigated from, both, the physical and 
mathematical viewpoints. For instance, the probability distribution of the time spent on a link 
(the so-called occupation time) was first studied by P Levy [13] who considered the time spent 
on an infinite half-line by a one-dimensional brownian motion stopped at some time t. This 
work allowed Levy to discover in 1939 one of his famous arc-sine laws [15]. Since that time, this 
result has been generalized to a star-graph [16] and also to a quite general graph [17]. Local 
time distributions have also been obtained in [IS]. 

It has been pointed out since a long time that first-passage times and, more generally, occupation 
times are of special interest in the context of reaction-diffusion processes |19[ [20] . Computations 
of such quantities in the presence of a constant external field have already been performed for 
one-dimensional systems with absorbing points (see, for example, [21J). This was done with the 
help of a linear Fokker-Planck equation [19] [22] . 

The purpose of the present work is to extend those results on a general graph with some absorbing 
vertices. We will consider a brownian particle diffusing with a spatially-dependent diffusion 
constant and subjected to a drift that is defined in every point of each link. The paper is 
organized as follows. In section 2, we present the notations that will be used throughout the 
paper. We discuss the boundary conditions to be used at each vertex in section 3 and, also, in 
the Appendices. More precisely, we analyse in details specific graphs in the Appendices A and 
B. The obtained results allow to deal with a general graph in Appendix C. Section 4 is devoted 
to the computation of the average time spent, before absorption, by a brownian particle on a 
part of the graph. In this section, we also calculate the Laplace transform of the joint law of the 
occupation times on each link. In the following section, we present additional results, especially 
concerning conditional and splitting probabilities. Various examples are discussed all along the 
different sections. Finally a brief summary is given in section 6. 

2 Definitions and notations 

Let us consider a general graph Q made of V vertices linked by B bonds of finite lengths. On 
each bond [a(3], of length l a p, we define the coordinate x a p that runs from (vertex a) to l a p 
(vertex (5). (We have, of course, xp a = l a p — x a p). 

Moreover, we suppose that, among all the vertices, N of them are absorbing. (A particle gets 
trapped if it reaches such a vertex). 

We will study the motion on Q of a brownian particle that starts at t = from some non- 
absorbing point x. The particle with a spatially-dependent diffusion constant D(x) is subjected 
to a drift U(x) defined on the bonds of Q. More precisely, D{x) and U(x) are differentiable 
functions of x on each link. In particular, on each link [a(3], the following limits Dr a p\ = 
^ m x a3 ->oD(x af 3), D' (a/3) = lim^^o di^ D ( x al3), U(af3) = lim Xa0 ^ o U(x a p) , are well defined. 
Such notations will be used extensively throughout the paper. 

The continuity properties of D(x) and U{x) at each vertex will be discussed in the following 
section. 

We also specify the motion of the particle when it reaches some vertex a. Let us call 
(i = 1, . . . , m a ) the nearest neighbours of a. We assume that the particle will come out towards 
Pi with some arbitrary probability p a ^ (SiPaft = 1 _ see [16] for a rigourous mathematical 
definition). Of course, = if A is an absorbing vertex or if [X/j] is not a bond of Q. 

Let P(yt\x0) be the probability density to find the particle at point y at time t (P(y0\x0) = 
5(y — x)). It satisfies on each link [af3] the backward and forward Fokker-Planck equations: 
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Figure 1: The vertex a with its nearest-neighbours /3j, i = 1, ...,m a ; each link is discretized with 
steps of lengths Ax. 
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L + (x a/3 )P(yt\x a f30) 



(1) 



dy, 



a/3 



dy a p 



dy a p 



L{y a p)P{y a (3t\xti) 



P{afi) wm mean lim Xc(/3 ^o P{yt\ x af$) when we use the backward Fokker-Planck equation and 
lim^^o P(y a/ 3t\x0) when we use the forward equation. The derivatives P'^p) wm ^ e defined 
in a similar way. 



3 Boundary conditions 



Let us define the following two situations that can occur at some non absorbing vertex a: 

(A) : D(x) is continuous but U(x) is not 

(B) : U(x) is continuous but D(x) is not 



The main purpose of this section is to establish the boundary conditions for P and its derivatives 
that result from those discontinuities. We will not consider the case when both U{x) and D{x) 
are discontinuous at the same vertex because, in our opinion, it is ill-defined. 

3.1 Backward Fokker-Planck equation 

Let us start by considering a graph Q with U(x) and D(x) constant (standard brownian motion). 
In Figure [H we display a given vertex a and its nearest-neighbours i = 1, ...,m a on Q. 

The transition probabilities from a, p a Pi, are not supposed to be all equal. In order to establish 
the boundary conditions for the backward equation, we discretize all the links [a/3i] (and also 
the time) with steps of length Ax (resp. At). It is easy to realize that: 
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Ilia 



P(y,(N + l)At\a,0) = £ p a0i P(y, NAt\ Xi , 0) (3) 



i=l 



P(y,(N+l)At\xi,0) = ^P(y,NAt\a,0) + ^P(y,NAt\x' i ,0) (4) 
Taking the limit Ax -> 0, At oc (Ax) 2 , AT -> oo, A^At = t, we obtain, with (gj): 



= ^2/, *l«> 0) + 2 P («ft) ( 5 ) 



Thus 

P{afk) =P(y,t\a,0) Vi (6) 

This shows that P is continuous in a. 
Moreover, expanding ([3]) at order Ax, we get: 

m a /rn a \ 

P(y, t|a,0) = $><*ft P (aft) + Ax $> aft P ( ' Q/3j) + 0((Ax) 2 ) (7) 
i=l \t=i / 

With © and ^SPaft = 1, we show that: 

nia 

E P U) = ( 8 ) 

i=l 



On the other hand, for = p Q/ g 2 = ... = l/m a , the equation ([I]) on the link [a/3] can be 
written in the form: 



^ ( 7 ^^e-*^)p(yt|x Q/3 0)') = f e -*M-^-P(yt|x^0) 



0) 



fan dU(x') dx' . , 

xo is some point on the graph. 

We are aware that $(x a p) could be multi-valued because, in general, a graph is multiply- 
connectecQ- However, this is not the case if we restrict ourselves to the vicinity of the given 
vertex a. Choosing xo located on a link starting from a, the integral in ()10p involves in a unique 
way, at most, two integrals along links starting from a. It is well defined if U{x) and D{x) 
are not discontinuous at the same point. So, with this definition of $(x Q/ g), the equation ([9]) is 
well-suited to study the boundary conditions at vertex a. 

Let us consider the case (A). 

We assume, first, that the p Q &'s are all equal. 



It can also occur that is not well defined, for instance, when U(x) and D(x) are both discontinuous 

1 dU(x') 

at the same point xi. Indeed, — is then proportional to 8{x — xi) / D(xi); thus, this quantity is not 

D(x) ox 

defined if D(x) is discontinuous at xi. 
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Integrating ([9]) in the vicinity of a, we get, with D{x) continuous at a: 

E e- U ^ /D(a) PU) = (11) 
»=i 

Moreover, P must be continuous at a if we want the quantity e _$ -r— to be properly defined. 

ox 

In the Appendix A, those boundary conditions are directly established, for the Laplace Transform 
of P, on a simple graph. 

Still for the case (A) , let us now assume that the p a 8 i 's are not all equal. 

In Appendix C, we establish the following boundary condition that must hold for a general 
graph: 

E PoA e- WW P^ ft) = (12) 
i=i 

We also show in this Appendix that P is continuous in a. 

Turning to the case (B), we can follow the same line as for (A) and use ([9]) and also the 
Appendices A and B. 

Remark that, when the PafcS (* = 1,2, ...,m a ) are all equal, the integration of (|9|) on the 
vicinity of a will not produce an exponential factor as in (|11|) . This is because, this time, 
U(x) is continuous in a. This fact is confirmed by the direct computations performed in the 
Appendices. Finally, with the Pafc's not all equal, we get: 

E P U) = (!3) 

i=l 

Moreover, P is continuous for the cases (A) and (B). 

In summary, for the backward Fokker-Planck equation, the condition on the derivatives can be 
written: 

E p'^ p U) = a*) 

i=l 

case (A) (15) 
case (B) (16) 

This notation will show especially useful in the following sections where the backward Fokker- 
Planck equation is widely used. 



with 
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3.2 Forward Fokker-Planck equation 



Coming back to Figure Q] (with Xi, x\ and Ax replaced by y^, y\ and Ay), we consider the 
discretized version of the forward equation with D{x) and U(x) constant and, also, the Waft's 
not all equal. We have, on the link [a/3i\ 

P( yi ,(N + l)At\x,0) = p a p i P(a,NAt\x,0) + ±P(y' i ,NAt\x,0) (17) 

P(a,(N + l)At\x,0) = P(yi,NAt\x,0) (18) 

z i=i 

With the limit Ay -> 0, At oc (Ay) 2 , N -> oo, iVAt = t, (JT7D leads to: 

i°(aA) = PafrP(u, t\x, 0) + -P( aft ) (19) 



Thus 



iW) _ P («/3 2 ) . F («/W _ 2 p(a,t|x,0) (20) 



We see that, in general, P is not continuous in a. 
Moreover, expanding (fTHl) at order Ay, we get: 



P(a,t\x,0) = ^ E ( + A ^ P U)) + °(( A ^) 2 ) ( 21 ) 

i=l 

With (|2U|) . we can write: 

E F U) = C 22 ) 

i=l 

So, the current conservation doesn't involve the p a /3i' s - 

Now, for Pafo = p a (3 2 = ■■■ = l/m a , the equation ([2]) on the link [a(3] can be written: 

^-P(y a0 t\xO) = ^(e^M^(D(y a ^Mp(y aP t\xO))) (23) 
ry a0 dU(v') dv' 



yo is some point on the graph; the discussion of section 13.11 for the definition of is, of course, 
still relevant. 

Following the same lines as in section [3~T| we get the current conservation at each non absorbing 
vertex a: 

Y,(( DP )' + U'P) {aM = (25) 
i=l 

if a is not the starting point. Otherwise, the right-hand side of ([25]) should be replaced by 
-6(t). 

(|25p doesn't depend on the continuity properties of D(x) and U(x). 
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Now, let us consider the case (A) and call D(a) the diffusion constant at a. 
When p a $ x = ... = p a p m , following [22], we can show that 

e P(aft) = eWP (aft) = ... = e ^5 F (a/3mJ (26) 
When the Pap^s are not all equal, we use the same approach as the one of Appendix C and get: 

(A) : e"°M" ( ll} = e "ow — = ... = e D <"> 1 ^ mtJ (27) 
Similarly, for the case (B), we obtain: 

(B) • D (a^)— — = D {a/ 3 2) —— = ... = £>(o/3 ma )— (28) 

Potpi Pa/3 2 Pa[3 ma 

When p a/ 3 1 = . . . = p a f3 ma > it is worth to notice that the conditions ([27]) and ([28]) are simply 
obtained if we want the quantity — (^D(y a p)e®( yal3 ^ P(y a pt\xO)^J , that appears in (j23j) . to be 
properly defined. 

In the following section, we will check the consistency of those boundary conditions (for the 
Backward and the Forward equations) on several examples. 

4 Occupation times 

4.1 Mean residence time 

Let us first study the average time, (r(x)), spent on a part V of Q by the particle before 
absorption. We have [21j : 

/•OO P 

(t(x)}= / dt / dyP(yt\xO) (29) 
Jo Jv 

{{t{x)) is the sum of the infinitesimal residence times weighted by the probability of presence 
in T> at time t). 

With ([I]), we get for (t(x)) the equation: 

L + (x a/3 ) {r{x a p)) = -l v (x a p) (30) 

lx>(x) is the characteristic function of the domain T>: 1t>(%) = 1 if x G P, = otherwise. 
The solution writes: 

(T(x a p)) = / du af3 <j>(u a p) + A a p / du a p I(u a/3 ) + B a/3 (31) 
Jo Jo 

with 

4>{x aj3 ) = - I{x a p) / \ -77 t- (32) 

Wo D(z a p) I(z a p) J 



I{x a p) = exp ^ 



<*0 dz a/3 <9C7 
D(z aj3 ) dz aj3 



(33) 
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The constants A a p and B a p are determined by imposing the boundary conditions at each vertex. 
Continuity implies: 

lim Mxap.)) = (T a ) = B a/3i (34) 
and lim (r(x Q/3 J) = (r ft ) = K^ p .) + A a/3i J (aft) + (35) 

where 

K (a/3 ) = / du Q/3 0(n Q/9 ) (36) 
o 

J (Qii g) = / du a/3 I{u a p) (37) 

J 

In general, K (a/3) ^ K (/3q ), J( Q/3) / J (/3q ). 

Of course, (rx) = if A is absorbing. So, we will determine (r a ) only for a non-absorbing. In 
such a vertex, the current conservation leads to ()14til6j) : 

£ P ' Ql3i A aPt = (38) 

i 

The equations for the (T a )'s can be written in a matrix form and finally the average time spent 
on D, before absorption, by a particle starting from a is given by: 

(r ) - (39) 
(Ta) "detM (<j9j 

M and Mi are two (V — N) x (V — N) matrices with the elements: 

M vi = J2t^ ( 4 °) 

Afy = -A- (41) 

J (ij) 

(In (|40|4ip , i and j run only over non-absorbing vertices but m labels all kinds of vertices) . 

th 

Mi = M except for the a column: 

m*=-Y,*JT& (42) 

We observe that, for a graph without any absorbing vertex, we have Mij = Vz. Thus, 
det M vanishes and (r a ) becomes infinite as expected. 

Remark also that, when D = Q, r becomes the survival time or, accordingly, the first-passage 
time in any absorbing vertex. 

4.1.1 Examples 

EXAMPLE 1. 

As a first example, let us consider the case T> = Q for the graphs of Figure [2j 
- in a), the vertex 1 is absorbing and the particle starts from 0. We get 

(r ) = - K {01) (43) 
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n 



Figure 2: Three simple graphs that are used to compute survival times (parts a) and b)) or 
covering time (c)). In a) and b), the vertex 1 is absorbing. The brownian particle starts from 
(parts a) and c)) and from A (part b)). For further explanations, see the text. 



AU 




X 0X X 01 

Figure 3: The potential is symmetric around except for a step of magnitude AU. 
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in b), 1 is still absorbing and the particle starts from A. Its probability transition in is 
Pqi = q (po\ = 1 — q). Moreover, we assume that Zqi = ^OAi U(xqi) = U(xq\) + AU (a step 
of magnitude AU occurs in 0, see Figure [3]) and D(xqi) = D(xq\). The mean survival time is 
given by: 



(n(?)> 



K, 



(01) 



(10) 



+ (1 



AU 

q) e D( ) 



(44) 



1 



c) represents a symmetric star with n legs of length I originating from (poi = — and 

n 

U(xoi) = U(xoj) (no step in this time) , D(xQi) = D{xQj) = 1, . . . ,n). Let us compute 
the average covering time (r c ) (smallest time necessary to reach all the points of the star at least 
once). With the notations of (pl3j) and (j4"4"|) . we get, with AU = 0: 



(45) 



n-l 



(rc) = (to) + £ - */«)> 



k=X 



n-l 



-%i) - (^(oi) + K{io)) 



k=l 



oo 



(46) 
(47) 



~ — (i^(oi) + -fiT(io)) Inn when n 
(Recall that for random walks on a star, (r c ) ~2n Inn [23J) 
EXAMPLE 2. 

Let us consider the graph of Figure [H Among the five vertices, two are absorbing (3 and 4). The 
starting point is 0. U(x) is discontinuous in vertices 1 and 2, D(x) is discontinuous in vertex 
0. All the links have the same length I. With the backward equation, we get the mean survival 
time (ESI): 



(T0) = 1 [ ^ + D- 2 



PmPlZ 



+ 



P02P24 



(48) 



vPi3 + Pioe^- u )l D ^ p 2 4 + P2oe( U2 - u ^ D2 . 

Let us now compute (tq) with the forward equation. Defining V(y) = Jq° dt P(yt\00), we have 

(49) 
(50) 
(51) 
(52) 



(53) 

(54) 
(55) 



on the link 


[31] : 


V{y zl ) = 


ax 2/31 + a-2 


on the link 


[10] : 


V(y x) = 


bx yox + h 


on the link 


[02] : 


V(y 02 ) = 


Cl 2/02 + c 2 


on the link 


[24] : 


V(y 42 ) = 


dx 2/42 + d 2 


with the boundary conditions 








P(3i) = o v m = 








Pox 


P02 






e U2/D 2V 


m _ ^ 


^(20) 


P13 P10 


P2A 


P20 



^(13) + ^(10) 



^(24) + ^(20) 



D 2 V' (t 



(02) 



(01) 



Solving those conditions and computing (to) = JgV(y)dy = (b 2 + c 2 )l, we recover the result 
equation (j48l) . 



EXAMPLE 3. 
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3 10 2 4 

G • • • O 

I I I I 



U(x) 





Figure 4: The graph Q consists in five vertices (3 and 4: absorbing, 0: starting point), U(x) is 
discontinuous in 1 and 2 and D{x) is discontinuous in 0. 

For the graph of Figure [5l the starting point is 0, a pure reflection occurs in 1 (pio = 1) and the 
vertex 2 is absorbing. The potential is U{xq2) = a £02> ^(^oi) = — a xq\. With the notations of 
the Figure and using the backward result (I39p . we get the mean residence time: 



on the link [10] (r ) = ( e ah/D2 - l) ( e all/Dl - l) (56) 

P02 a 2 v ' v ' 

on the link [02] (r > = ^ ( e ah ' D2 - I - ^f] (57) 

a 2 \ D 2 J 

Taking the limit a — > (no drift), we obtain: 

on the link [10] (r ) = — ^ (58) 

P02 Di 

I 2 

on the link [02] (r ) = (59) 

Turning to the forward equation, still with V(y) = dt P(y£|00), we have 

on the link [10] : P(m) = cy e ayoi/Dl + c 2 (60) 
on the link [02] : P(y 02 ) = c 3 e" ayo2/D2 + c 4 (61) 

with the conditions 



^ 2 o) = o ^m = ^m Mlo) + ffl p (10) = o (62) 

v ; P01 P02 K ' v ; 

D 2 V' m + £>i7f 01 ) + a (^(02) " P(oi)) = -1 (63) 



Computing (t > = Jj,P(y) dy = ci^ (e a/l / Dl - l) +c 2 /i for the link [10] , = cs^ 2 - (l - e" ' 2 ^ 2 ) + 
C4/2 for the link [02], we recover the solution (|56l57p . 
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Figure 5: The graph Q consists in three vertices (2: absorbing, 0: starting point, pure reflection 
in 1); the drift U(x) is linear and D(x) is discontinuous in 0. 

4.2 Occupation times distribution 

Still with the same conditions (drift, spatially-dependent diffusion constant, absorbing ver- 
tices,. . . ), we now study a Brownian motion stopped at t and call T Q( g the time spent up to 
t on the link [ot(3]. 

T a p is, of course, a random variable depending on t and we will call Vt{{T a p}) the joint distri- 
bution of the T a ps. 

In the following, we will focus our attention on the quantity 

S(t\x0) = [Y[dT a pV t ({T a p})e-^^ (64) 

[a/3] 

-£rc]*« r «(a;)\ (65) 



where the &j's are positive constants and (. . . (x)) stands for averaging over all the brownian 
trajectories starting from x and developing up to t in the presence of the drift. 

In order to stick to S(t\x0), we slightly modify the backward equation ([1]) by adding, on each 
link [a/3], a loss term proportional to £ a( g may be interpreted as a reaction rate per unit 
length and time). Thus, we consider the equation: 

dO 

-Tj^ = (L + {x a /3) - £ a p) Q (66) 
where Q(yt\x0) represents the probability of finding the particle in y at time t, each path from 

— ^ t 

x to y being, this time, weig hted by a factor e ^W]^ 1 ^ . (Notice that Q(yt\x0) = 5(y - x) 
if x is absorbing). 
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It is now easy to realize that we have the relationship: 

S(t\xO) = f dy Q(yt\xO) 

Jg 

In the following, we will be especially interested in the Laplace Transform of S(t\xO): 

/■oo 

S(7|x0) = S(x) = / dt e" 7 * S(t\xO) 



On the bond [a[3], S(x a p) satisfies the following equation (7^ = 7 + 

(L + (x a p) - 7 Wj9 ) S = -1 

Setting S{x a 3) = \- Q(x a n) and performing the transformation 

la/3 

(I is defined in (|33p ). we are left with the following equation for \'- 



(67) 



(68) 



(69) 



dx l(3 



1 dU{x. 



a/3) 



4D \ dx, 



D 



d I 1 dU{x a f 



a/3 



dx a/ 3 \2D dx 



a/3 



+ la(3 



(70) 



X = (71) 



Let us call Xa/3 and Xf3a two solutions of (fTTj) such that Xap( a ) = Xfia(P) = L Xfiaip) 
Xa/3(P) = 0. So, 5 writes: 



S = — !— + A a R exp 



+ A^q, exp / 



x <*f> dz af5 dU{z a p) 
2D(z a p) dz aj3 

<*0 dz afi dU(z a/3 ) 
af3 2D(z Q/3 ) dz a /3 



Xa/3 (X a/?) 
X/3a{x a /3) 



(72) 



The constants Aij are determined by imposing the boundary conditions. Continuity of S at 
each vertex a implies: 



lim n S(x a p.) = S a = — !— 



A 



a/3, 



Moreover, if a is absorbing then = V [Xfi] and, from (164j) and (168p . we get: 



S -I ' 



A. 



7 Taft 



a/3, 



(73) 



(74) 



Following the same lines as for {r a ), it is easy to show that the S Q 7 s (a non-absorbing) satisfy 
a matrix equation with the solution: 



S a 



det Ri 
deti? 



i? and i?i are again (V — N) x (V — iV) matrices with elements 

Rii = Pirrfiim 



Ri 



PijWij 



(75) 

(76) 
(77) 
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the quantities Cj m and Wij being given by: 



' a dU 



%im~ *o 1 dxi m 2Z)(a;j m ) J 



9C/ 



^ = ^i£4feH*«"^f (79) 



th 

R\ = R except for the or column: 

( R i)i« = J2^-( C v+ W ij)-l E PikWik (80) 
i llJ 7 fcabs. 

(The last summation is performed only over absorbing vertices). 

1 

Setting £ij = 0, we establish that det R\ = — det R, showing that the probability distribution 

7 

T > t({T a /3}) is properly normalized (see (f64"|) and ([68]) ). 

Now, let us call T' a a the time spent on [a/5] before absorption. Thus, for a particle starting from 
a, we can write 

r; -E W ]^( a A = ^/ e ~£w]* yTy (a)\ (81) 
and from (|65l68p we deduce 

e" ^W] ^ (a)\ = lim ( 7 5 Q ) (82) 
Finally, we get: 

with it!^ = lim R and i?^ = R^ except for the column that is given by 



7->0 



HU) holds because 



(4 0) ). =-hm J2 PikWik (84) 
k abs. 



This is obvious when 7^ 0. But when = 0, this is still true because, in that case, 

lim (C i:j + ) = 0. 
7 — >o 

Remark that (|84p implies that ^ T «(a)y = if £/ has no absorbing vertex (in that 

case, T(j = 00). 

Setting £ij = £ and ^ T/^ = t (survival time), ([83]) gives the expression of the Laplace Transform 

[ij] 

of the probability distribution of r. 
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4.2.1 Example 



e ~to\ = _ IMU (85) 



To illustrate this work, let us go back to the example of Figure [2] and study the survival time 
(parts a) and b) - with AU = 0) and covering time (part c)) distributions. 

With a), we get 

Coi 

(W a p, C a p,. . . are computed with j a p = £). 
b) leads to 

(e-^(g))= - ^ •', , (SO) 

with Z 



Z-l+q 

CoiCio 



For the covering time r c in c) we get (with the notations of (1851) and ([86 



n-1 

e- fTc ) = - (e~ 5T0 ) J] (e-^{l-k/n)) (87) 

fe=l 



Wb! (n-1)! n fr 1 f88l 

k=l n 

Computation of the first moments shows that when n — > oo: 

(t c ) oc n Inn in agreement with (|47p (89) 
(r c 2 )-(r c ) 2 oc n 2 (90) 

Then the probability distribution of the scaling variable r c j (r c ) becomes more and more peaked 
at its mean value when n — > oo: 

P = ^„^oo <5(X - 1) (91) 



5 First passage times 

This section is essentially devoted to a detailed study of exit (or survival) times. First, we 
will compute exit times distributions and, after, in the last parts, we will focus on quantities 
concerning exit through a given absorbing vertex. 

5.1 Exit times distribution 

The density of exit time V(t\x0) (without specifying the absorbing vertex) is defined through 
the loss of probability: 

V(t\xO)dt = [ dyP(yt\xO)-[ dy P(y, t + dt|x0) (92) 
Jg Jg 

Thus: V(t\xO) = ~ [ dy P(yt\x0) = T V^t\x0) (93) 

at Jg , 

fi abs. 
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V^(t\xO) is the density of exit time by the absorbing vertex fj,. 
With the forward Fokker-Planck equation, (j93[) leads to: 



V(t\xO) = -J dy L(y)P(yt\xO) = £ (94) 

/j,abs.,i 

( d dU \ 
= lim l-£-(DP) + - P) (95) 

lim (d^-) (96) 



because PupA = when \i is absorbing. (The /%'s are the nearest-neighbours of fi on Q.) 
We deduce that Vn(t\xO) is the probability current at vertex fj,: V^^xO) = -^tft- 

i 

Let us compute the Laplace transform, V^xO), of V(t\xO). With ([93]) . we get: 

P( 7 |xO) = J^Me-V^-lLj^ dyP(yt\xO)^ = £ V^\xO) (97) 



; abs. 

= 1 - 7 / die" 7 * / dy P(yibO) (98) 

JO JQ 

The backward Fokker-Planck equation gives: 

L+V^xO) = 7^(7|x0) (99) 
and L + V^\xO) = jV^\xO) (100) 

Defining V^ a = lim 7 3 ^( 7 |x a i0), we remark that 

for an absorbing vertex A : Pp.,\ = <^A (101) 

This is because ^(t|A0) = bit) if A = fj, (102) 

= otherwise (103) 

Following section H~2j we readily find, for a particle starting from a (non absorbing), the Laplace 
transform of the exit (by the vertex /i) time distribution: 

v -° = -35*- (104) 

where R^'°^ = R except for the cfi 1 column: 

? (jj,,a) 



R * )ia = (105) 

R and are defined in (|76M79p and computed, here, with = 0. (p^ is defined in (|15H16p ). 
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5.2 Splitting probabilities 

For a particle starting at x, we define the splitting probability ILj(x) as the probability of ever 
reaching the absorbing vertex \x (rather than any other absorbing vertex). Such a quantity 
has already been considered for one-dimensional systems [23] (see also |25j for extensions to 
higher-dimensional systems). 

We have, obviously: 

n„(x) = / dt V„(t\xO) = = 0|x0) (106) 
Jo 

Setting 7 = in (jlOOp we see that II M (x) satisfies, on [a/3], the backward equation: 

L + {x a p)H^x aP ) = (107) 

The probability, ILj jQ , for a particle starting from the vertex a to be absorbed by fi, is defined 
as U^a = lim U^Xai). 

With (1106j) and (llOip . it is easy to realize that 

for an absorbing vertex A : IL^a = (108) 

Following the same lines as previously (see section 14. ip , we find that ILj ja (a non absorbing) is 
again written as the ratio of two determinants: 



detM 



= ~TrL ( io9 ) 

(M>°0 _ A/ r + u Q ^,th 



detM 

where M^*'" 1 ' = M except for the a lli column: 



(M^ a) ) = (no) 



(M is defined in (1401411) and in (|57|) ). 

With simple manipulations on determinants, we check the normalization condition ^ ^ IT, a 

i abs. 

1. 



Let us, for one moment, comment the case when there is no drift (U(x) constant but D(x) 
variable, eventually discontinuous at some vertices). In that case p'^ = pij and, also, = kj. 
We conclude that the splitting probabilities don't depend on the varying diffusion constant 
when there is no drift. This fact can be understood in the following way. Let us consider a 
discretization of each link and a continuous time. Modifying the diffusion constant amounts to 
change the waiting time at each site of the discretized graph. But, this would not change the 
trajectories if there is no drift. Only the time spent is changed. Finally, the splitting probabilities 
remain unaffectecH by a change of D(x). 

2 For a graph Q without drift, we could expect, with the same argument, that the average time spent on a part 
D of Q would not depend on the diffusion constant on the rest, Q\D , of the graph. This is exactly what can be 
checked with the formulae (|39H42p of section fOl and, also, with the formulae (|58l59p of the example 3 of section 

mil 
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5.2.1 Example 



Let us consider a star-graph without drift. The root has mo neighbours, all absorbing. The 
ijiq links have lengths loi, i = 1, ...,mo. With (|109p . we obtain 




5.3 Conditional mean first passage times 

We now turn to the study of the conditional mean first passage time (^(x)}, which is defined 
as the mean exit time, given that exit is through the absorbing vertex [i (rather than any other 
absorbing vertex). We set {t^ a ) = lim (tfj,(x a i)}. 

Actually, it is simpler to first compute the quantity 9^{x) = n^(x) (i^(x)). 
Indeed, we have: 



X' 



e^x) = / dttV^(t\xO) (112) 
Jo 

rco Q poo 

L+B^x) = / dt t — r^t\x0) = - / dtV^xO) (113) 
JO ot Jo 

So, we get: L+O^x) = - U^x) (114) 
Moreover, for any absorbing vertex A, we get: = ILj^ {t^,x) = <5/x,a {tfi,n) = because 

Thus, comparing this equation with eq. (|30p . we find for a particle starting from the vertex a 

n - <w> = (us) 

where M^' a ^ = M except for the oft 1 column: 

(M^). = E^T^ (H6) 

V Ua m J (im) 



with 



s$ = / o "d„,^-) /;« d%i i|L ( U7) 



In this last equation, Il^(zij) has to be computed by equation 



IL^Zij) = n„,j + UflJ T U ^ I du„ I i n,, l (118) 



J(ij) Jo 
with IL^j given by eq. (|109p . 
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5.3.1 Example 



With the same star-graph as in section 15.2,11 (no drift) and a diffusion constant equal to Doi on 
each link [Oi], we get (|111I115P 

6 Summary 

We have used the Fokker-Planck Equation (mainly the backward one) to compute some quan- 
tities relevant for the study of a Brownian particle moving on a graph. This particle, moving 
with a varying diffusion constant, is also subjected to a drift. 

Being aware that this point is somewhat controversial (see, for instance, [26])) we have discussed 
in great details the boundary conditions in vertices where either U(x) or D{x) is discontinuous. 
Two Appendices support the study of those boundary conditions. 

Those preliminaries allowed us to deal with various quantities such as the mean residence time, 
the occupation times probabilities, the splitting probabilities and the conditional mean first 
passage time. In each case, we have established general formulae that can be used for all kinds 
of graphs. 

Finally, we have checked, on several examples, the consistency of the boundary conditions we 
have put forward. 



A Backward Equation: boundary conditions when poi = P02 

For a direct computation of the boundary conditions, it is simpler to consider the Laplace 
Transform: 

S(j\xO) = / dt e" 7 * / dyP(yt\xO) (120) 
jo Jg 

that satisfies the backward equation (5(71x0) = S(x)): 

(L+- 7 )S(x) = -1 (121) 

First, we want to study the case (A) 

For the simple graph of Figure [6] a), D(x) is continuous and U(x) is discontinuous at the vertex 
0. In particular f(oi) = U\ ^ £^(02) = U2',D(0) = D. Moreover, vertices 1 and 2 are absorbing 
and, in 0, the transition probabilities are equal: poi = P02 = 1/2. 

Now, let us add to this graph, the link [03] of length L (see Figure [6] b)) in such a way that 
nothing is changed for the rest (for example, the potential on the link [02] of the original 
graph is the same as the one on the link [32] of the modified graph, ...). On [03], we define 
^(^03) = (~ 2 tt X j x 03 + U\ and Z?(xo3) = D. When L — ► 0, we recover Figure [6] a). 

Working with this modified graph, the solution of (I121D writes: 
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Figure 6: In part a) a simple graph with the absorbing vertices 1 and 2; the potential is 
discontinuous in 0. In part b), we add the link [03] (length L) and, on this link, a potential that 
interpolates linearly between U\ and U 2 \ D{x) is constant on [03]. 



on the link [10] S(x\q) = I/7 + ai^i(xio) + a 2 (f> 2 (x w ) (122) 
on the link [03] S(x 03 ) = I/7 + Cl e r+X ° 3 + c 2 e r - x ° 3 (123) 
on the link [32] S(x 32 ) = I/7 + hipifaz) + b 2 ^ 2 (x 32 ) ( 124 ) 



with 



A±y/A> + 4ryD 
r± = w (125) 

A = (126) 
L 

0i 5 2 and -01,2 are solutions of the equation (L + — j)(f> = that satisfy (4>( a /3) = li m x a ^->o 4>( x ap))'- 

4>\ (10) =° <£i (01) = 1 ^2 (10) = 1 ^2 (01) = ( 127 ) 

V>1 (32) = 1 V>1 (23) =0 ^2 (32) =0 ^2 (23) = 1 ( 12 8) 

U{x) and D(x) are continuous everywhere on the graph. So, we write standard boundary 
conditions at the vertices (continuity of S and its derivative at vertices and 3; S = at 
vertices 1 and 2). We obtain (<j>[ a/3 ) = lim^^o 

a 2 = b 2 = -- ai = Cl + c 2 h = Cl e r + L + c 2 e r - L (129) 
7 
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£(32) 


i + 6l 


-L^O 1 


^(01) 




13 (32) 


6 lV4 (32) 


+ & 2 V4 (32) 


or 

*>1) 


(01) 


+ a 2 <t>' 2 (01) 



a l^'l (01) + a 2^2 (01) + C l r + + C 2 r - = ( 13 °) 

fel^i (32) + (32) " c ir+ e r + L - c 2 r.e r - L = (131) 
After some algebra, we get the results: 

(132) 

>L ^ _ e {U*-Ui)/D (133) 

When L — » 0, the vertex 3 moves to and equations (|132p and (| 133j) give, for the original graph: 

£(oi) = £(02) (134) 
$' m e~ Vm,m + S{ 02) e- u ^/ D ^=0 (135) 

In particular, we observe that S is continuous at and that exponential factors appear in the 
condition involving the derivatives. 

Let us now turn to the case (B), still with the same graph and poi = P02 = 1/2. 

This time (Figure [7] a)), D(x) is discontinuous at and U(x) is continuous. We set: -D(oi) = 
D x + D m = D 2 , U(0) = U. 

As we already did for the case (A), we modify the graph and obtain Figure [7] b). Between vertices 

\ 2 



and 3, we choose: U(xq^) = U and Z?(x 3) = — ^"~7j J x ° 3 V^iJ = ( ax 03 + b) 2 . 

Thus, in the modified graph, D{x) and U(x) are continuous everywhere on the graph. 

On the links [10] and [32], the solution of equation (fT2"Tj) is still given by (fT22|) and ([HID ( with 
the conditions (fT2Tj) and (fT28j) ). But, on the link [03], (H23J) has to be replaced by 

S(x 03 ) = 1/7 + ci (a x 03 + 6) A + + c 2 (a x 03 + &) A ~ (136) 
with: A± = ±±±^1 + ^1 (137) 



L 

Standard boundary conditions imply: 



(138) 



a 2 = b 2 = -- a x = Cl D X+/2 + c 2 D\- 12 b x = Cl D X+/2 + c 2 D^ /2 (139) 
7 

aicP[ m +a 2 ^ 2m +a[ Cl X + D[ X+ - 1)/2 + c 2 X^d[ X -- 1)/2 ) = (140) 

bM^+M^s^-^^lA+^^ + QA^-- 1 ^ 2 ) = (141) 



Finally, we get: 



?(32) " " + h >^o 1 (142) 



£(01) 7 + °1 
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a) b) 



1 2 1 3 2 




-l 1 1 2 -I j L L+l 2 



Figure 7: a) The same graph as in Figure [6] a) but, this time, D{x) is discontinuous in 0; in 
part b), a parabolic interpolation is used to make D(x) continuous; U{x) is constant on the link 
[03]. 



5 (32) Mi (32) + M>2 (32) . „ ... 
S (01) ai<f) l (01) + a2 ^2 (01) 

When L — * 0, we get, for the original graph, S continuous in and 

3(oi) + s'm = ( 144 ) 



B Backward Equation: boundary conditions when poi P02 

We want to study the case (A) (U(x) discontinuous in 0) for the graph of Figure [8] a) with, this 
time, pox / p 02 . 

Modifying this graph, we obtain a new one, consisting in five vertices, displayed in Figure [8] 
b). Vertices 1 and 2 are still absorbing. We set = p^o = p^o = p%2 = 1/2 but po4 P03 ( 
Poa = Poi (original graph) and p 3 = P02 (original graph)). 

U(x)(resp. D{x)) is set equal to some constant C/(resp. D) on the added links [40] and [03]. 
U (x) is discontinuous at vertices 4 and 3. With the modified graph, we can take advantage of 
the computations of the boundary conditions performed at the beginning of section [3] and also 
in Appendix A. 

The solution of equation (|12ip writes: 
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a) 







b) 



4 3 

lj L V l 2 




-L-lj 



U(x). 

u 



-L V 



L'+l, 



-h 



D(x) 





l 7 -L-L 



-L U 



L'+l 



Figure 8: a) The same graph as in Figure [6] a) but, this time, poi 7^ P02] in b), U(x) and D(x) 
are constant on the added links [40] and [03]; D(x) is still continuous on the modified graph. 



on the link [14] S(xu) = 1/7 + ai</>i(xi4) + a24> 2 {xiA) (145) 



on the link [40] S(x 40 ) = 1/7 + h S mh(^j/Dx 40 ) + b 2 smh(^/D(L - x 40 )) (146) 
on the link [03] S(x 03 ) = I/7 + a sivh(^/Dx 03 ) + c 2 a\ria(ypy/D{L' - x 03 ))(147) 
on the link [32] S(x 32 ) = 1/7 + di ^1(^32 )+ d 2 ^2(^32) (148) 
As before, ^2 and -01,2 are solutions of the equation (L + — 7)^ = that satisfy: 

<t>\ (14) =0 4>i (41) = 1 </>2 (14) = 1 <t>2 (41) = (149) 
^1 (32) = 1 V>1 (23) =0 ^2 (32) =0 1p 2 (23) = 1 (150) 

The boundary conditions in are given by the beginning of section 13.11 S is continuous and 

£>04 S(04) + P03 5(03) = (!51) 

For vertices 4 and 3, we use the result of Appendix A (case (A)). S is again continuous and 

S' m e~^ D + S' m e- U ' D = (152) 
S' m e~ U/D + S' m e~ u */ D = (153) 

We get the relationships: 
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1 



a 2 = d 2 = -- ai = b 2 smh.(^j/DL) d% = c x sinh(y j/DL') (154) 



h sinh(J~//DL) = c 2 sinhUj/DL') (155) 



P04 coBhty-y /DL) + 6 2 J + poa (^1 - c 2 cosh(^/DL')j = (156) 
(ai^i {41) + a 2 0' 2 (41) ) e"^ + yfyJB L - b 2 cosh(^/DL)\ e~ u / D = (157) 



(di^ 1{32) +d 2 ^' 2{32) )e- u ^ D - y^/D ( Cl cosM^tTdL') - c 2 ) = (158) 

Solving and taking the limit L, L' — ► 0, we are left with: 



•5(32) 7 + ^1 
5(41) 7 + °1 



5; 



(32) _ (32) + d 2^ 2 ( 32 ) U 2 -Ui 



e^(-M) (160) 

^(41) (41) + ^2 (41) V ^03 / 

When L, L' — > 0, the vertices 3 and 4 move to and we gelH: 

5(02) = 5( 01 ) (161) 
P02 5( 02) e " W<°> + Pol £ e -W>«9 = (162) 



We also observe that: 



5(03) 

~ >L,L'^0 

•3(30) 
0/ 

D (Q3) 

"& %L'^0 

O (30) 

Those relationships will show useful in Appendix C. 



(159) 



(163) 
(164) 



The solution for the case (B) (Figure 9 a), poi P02) is immediate. Indeed, we see that for 
the modified graph (Figure [9] b)), we have, formally, the same solution as equations (1145p - (|148|) 
and the same boundary conditions but, this time, with U\ = U 2 = U. It amounts to drop 
the exponential factors in (|152() and in all the equations that follow. Finally, we get that S is 
continuous and 

P02 5( 2) + P01 5( 01) = (165) 

Moreover, ([163]) and (fTM]) still hold. 

3 Recall that po4 = P01 (original graph) and po3 = P02 (original graph). 
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a) b) 



1 21 403 2 




-i 1 l 2 -L-lj -L L' L'+l 2 



Figure 9: a) The same graph as in Figure [7] a) but, this time, poi 7^ P02] in b), U(x) and D(x) 
are constant on the added links [40] and [03]; U(x) is still continuous on the modified graph. 

C Backward Equation: general boundary conditions 

For the case (A), let us assume that the p Q p^s are not all equal. In Figure [10] a) , where a given 
vertex a and its m a nearest-neighbours are shown, we suppose that U(x) is discontinuous in a 
(D(x) is continuous, D(x) = D(a) in a) . 

In part b), we slightly modify the graph along the same lines as in Appendix B. We add vertices 
[ii in such a way that each new link [/ij/3j] is identical to the original link [a(3i] (same length, same 
potential and diffusion constant). Moreover, in the added subgraph (Figure [10] b), heavy lines 
of lengths L) the potential and the diffusion constant are assumed to be constant (respectively 
equal to some value U and to D(a)). So, in the vicinity of a, the discontinuities of U(x) will 
occur at the /ij's (D(x) will be continuous in the same domain). Of course, for the transition 
probabilities from a, we choose p a ^ = PaPi- By taking the limit L — > 0, we will recover the 
original graph. 

Now, for the small subgraph where U(x) and D{x) are constant, we can take advantage of the 
result, equation ([8]), to write: 

E Paw P[ am) = (166) 

i=l 

Moreover, for the vertex fii, where p^ a = p ll ^ l = 1/2, we can use (fTTj) and also f|135[) (directly 
established in Appendix A) to get: 

<r U/D{a) PU«) + z- UM,D{a) PU) = (167) 
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a) 



b) 





Figure 10: a) the vertex a with its nearest-neighbours f3{, i = 1, ...,m a ; b) we have added the 
heavy lines where the potential is set equal to some constant U and the diffusion constant is 
equal to D(a); for the rest of the graph, nothing has been changed. For further explanations, 
see the text. 



Weighting with p a p.(= p atH ) and summing over i, we deduce: 

V T) a e- U (^i)/ D ^P' \ - - e -U/D{a) V- pi 



i=l 



\i=l 



Now, taking the limit L — * 0, we have from (|164p (Appendix B): 



pi 
pi 



y L^0 



(168) 



(169) 



In this limit, the vertex Hi moves to a and we recover the original graph. Finally, with 
(|166|168ll69p . we obtain, for the case (A), the boundary conditioro 



(170) 



Moreover, for the modified graph, P is continuous in a and in /Uj. From Appendix B, we know 
that (fT63j) : 

(171) 



(Mi a) 



P, 



1 



That's enough to conclude that, for the original graph, P is continuous in a. 
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